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ABSTRACT 

Weak lensing is commonly measured using shear through galaxy ellipticities, or using the effect of magni- 
fication bias on galaxy number densities. Here, we report on the first detection of weak lensing magnification 
with a new, independent technique using the distribution of galaxy sizes and magnitudes. These data come for 
free in galaxy surveys designed for measuring shear We present the magnification estimator and apply it to 
an X-ray selected sample of galaxy groups in the COSMOS HST survey. The measurement of the projected 
surface density I](r) is consistent with the shear measurements within the uncertainties, and has roughly 40% 
of the signal-to-noise of the latter We discuss systematic issues and challenges to reaUzing the potential of this 
new probe of weak lensing. 

Subject headings: gravitational lensing: weak — large-scale structure of universe — dark matter 
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1. INTRODUCTION 

Weak lensing has emerged as a powerful tool in cosmol- 
ogy and astrophysics (e.g., Refregier (2003b); Munshi et al. 
(2008); Hoekstra & Jain (2008); Massey et al. (2010); Schrab- 
back et al. (2010)). It probes the underlying total matter dis- 
tribution, and can be used to tie the observed distribution 
of galaxies to that of the dark matter, the dominant com- 
ponent of the cosmic matter budget. So far, weak lensing 
has primarily been measured using the shear, through the 
statistical correlation of galaxy ellipticities (see Bartelmann 
& Schneider (2001) for a review). In addition, weak lens- 
ing magnification has been detected through its effect on the 
number density of a flux-limited sample (magnification bias; 
Scranton et al. (2005); Hildebrandt et al. (2009)). Shear 7 
and convergence n, which are related to the magnification /i 
through jjL = [(1 -K-)+ |7|~]~', so that yU = 1 +2k to linear or- 
der, measure different properties of the density field. For an 
azimuthally symmetric lens. 
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where I](r) is the projected surface mass density of the lens, 
r is the physical transverse distance on the lens plane, and 
Ds, Di, Dis denote the angular diameter distance out to the 
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source, lens, and between lens and source, respectively. Note 
that the shear is proportional to the differential surface mass 
density AE, while the convergence probes I](r) itself. The 
different radial dependences of k, and 7 can be used to break 
degeneracies between the lens mass and density profile shape 
(Rozo & Schmidt 2010). 

The convergence k, has previously been measured using the 
observed number density of background sources. In a flux- 
limited survey, the observed number density is given by 



"obs = n[\ + 5„ + {5s-2)K\, 



(4) 



where fi is the mean number density, 5g denotes the intrin- 
sic fluctuation in galaxy density, and s = dlnii / dm^it for a 
sharp magnitude cut njcut- This can be generalized to surveys 
that are also limited by a size cut (Schmidt et al. 2009b). By 
cross-correlating foreground lenses with background sources 
widely separated in redshift, one can isolate this effect to mea- 
sure K. Note that the intrinsic fluctuations Sg act as noise, and 
any residual physical correlation with the lenses (e.g., due to 
photo-z uncertainties) will contaminate the signal. 

In this paper, we report on a detection of magnification 
around X-ray selected galaxy groups in COSMOS using 
galaxy sizes and fluxes. Related approaches have been pro- 
posed in Jain (2002); Bartelmann & Narayan (1995); Berlin 
& Lombardi (2006); see also Menard et al. (2010) who used 
individual quasar magnitudes, and Huff & Graves (201 1). We 
rely on the observed properties of individual galaxies and do 
not use the number density of galaxies as done in magnifica- 
tion bias measurements. One key advantage of this technique 
is that the contamination due to residual physical correlation 
with the lenses is likely much smaller than in magnification 
bias measurements, due to the weaker correlation of galaxy 
sizes and luminosities with environment (Maltby et al. (2010); 
Nair et al. (2010); Croton et al. (2005)). 

Throughout, we adopt a flat ACDM cosmology with h = 
0.72, il„, = 0.258 to calculate distances. Halo masses are de- 
fined through a mean interior density of 2QQp„h{z = 0.38), 
where z = 0.38 is the mean redshift of the lens sample. 

2. MAGNIFICATION ESTIMATOR 

The data set we use consists of a set of values {ii,-,m,,z,} for 
each galaxy, denoting the size, magnitude, and photometric 
redshift. Let us consider how these observables change under 
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lensing magnification. Tliroughout, we will work in the weak 
lensing regime. We can then write 



d = {l+riK)do 
m=mo + qK 

Z=Z0. 



(5) 
(6) 
(7) 



Here, a subscript denotes unlensed quantities, and in the 
last line we have assumed that photometric redshifts are not 
affected by lensing. In absence of noise and instrumental 
effects, the efficiency factors are i] = I and = -5/lnlO ~ 
-2.17, respectively (e.g., Schmidt et al. (2009b)). In prac- 
tice, these values can be different due to the instrument PSF, 
and other effects and in general depend on the properties of a 
given galaxy. 

In order to estimate k using Eqs. (5)-(6), we need to know 
the unlensed quantities do^niQ. These have some intrinsic dis- 
tribution, and given the known distribution we can infer n sta- 
tistically. Our starting point for the estimator is thus the joint 
distribution p^)(do,mo\z) of unlensed galaxy sizes and magni- 
tudes at fixed redshift; the estimator is designed to not use the 
observed number density of galaxies. We approximate po as 
Gaussian in m and x = Ind. Furthermore, we approximate the 
error distributions in m and x as Gaussian. This greatly sim- 
plifies the evaluation of the estimator, at the price of making 
it sub-optimal. The log-likelihood is then given by 
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The total variance in x and m is given by the sum in quadra- 
ture of the measurement uncertainty, intrinsic dispersion, and 
propagated photometric redshift uncertainty. CTx.meas, tT„, ,neas 
are measured as a function of d and m, respectively, through 
multiple images of the same galaxy from overlap regions of 
the HST ACS data (Leauthaud et al. 2007). The propagated 
redshift uncertainty is approximated as {a = x, m) 
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where z±,i are the upper and lower 68% confidence level val- 
ues for the redshift of galaxy /. The corrections e^, e„, to fj 
and q take into account the fact that a positive magnification 
pushes faint, small galaxies over the flux and size thresholds 
(lensing bias). A similar effect exists for shear, but only at 
second order (Schmidt et al. 2009a), while the magnification 
estimator is affected at leading order The correlation between 
intrinsic galaxy sizes and magnitudes is quantified by the cor- 
relation coefficient p. 

The maximum-likelihood estimate E and its estimated vari- 



ance are then given by 
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where 



Ai = - — \— I ^ (xi - x{zi)) + -f - (nii - m{zi)) 



P(Zi) 



[fiiinii - m{zi)) + qiixi - x(z,))] 
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The intrinsic dispersions o-.v.int, C/n.int as well as p are deter- 
mined by fitting a bivariate Gaussian distribution to all galax- 
ies before cuts on x, m in a redshift slice. We subtract the 
mean measurement error in x, m in quadrature from the total 
measured dispersion to obtain cr.i.int, dmM- 

The quantities x, m are determined by solving the equation 
E = applied to all galaxies in a given redshift slice. Since 
F w const within a narrow redshift slice, this yields 
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y^-, \ -^imi-m(zi))- ^^^'^ fiiinii-mizi)) ) =0. 

; 1 - P iZi) y ct;;,j crx,icrm,i J 

(19) 

X, m do not correspond to the true mean intrinsic quantities, 
since only galaxies passing cuts are used to measure them. 
We can then write 



a(z) = \ 
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(20) 
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where the sum runs over all galaxies passing cuts in the red- 
shift slice. 

In order to determine E j , £»» we repeat the measurement of 
X, m with magnitude and size cuts varied by ±Am, ±Ax, re- 
spectively, where Am = qinicudKo and Ax = ri{dcat)Ko, and we 
take kq = 0.02 (the value of e, does not depend significantly 
on kq). Then, 



£x(z) = 

£m(z) = 



x+{z)-x-{z) 

2kq 
m+(z)-m-{z) 

2ko 



(23) 
(24) 



where a± indicate quantities measured when using cuts fWcut ± 
Am and Xcut ± Ajc. In order to take into account interdepen- 
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Fig. 1 . — Upper panel: Average source galaxy log size x{z) as a function 
of redshift, measured using Eq. (20). The shaded band shows (Tvjnt- Lower 
panel: Same, for source galaxy magnitudes (m(z), cfm.int)- 

dencies, we iterate the measurement of x(z), m{z) until con- 
vergence. 

3. DATA 

Our data set consists of the HST ACS imaging of the COS- 
MOS field (Koekemoer et al. 2007), from which we take 
the galaxy sizes d (in pixels of 0.03") and magnitudes m in 
the F814W band. For the size d we use the variable-width 
Gaussian-filtered second moment of Rhodes et al. (2000), 
while magnitudes m are MAG_AUTO from SExtractor (Bertin 
& Arnouts 1996). We leave more sophisticated, optimized 
size and flux estimators for future work. The details of the re- 
duction pipeline are described in Leauthaud et al. (2007). The 
size estimator is PSF-deconvolved, and we thus expect sys- 
tematic effects of PSF variations to be negligible. In order to 
avoid any contamination by a background-dependence of the 
estimated sizes, we only use sky regions within the lower peak 
(mean pixel background < 0.0062) of the essentially bimodal 
background distribution in the COSMOS ACS data, compris- 
ing - 75% of the field (see Fig. 8 in Leauthaud et al. (2007)). 

For our source galaxy sample, we exclude galaxies in 
masked regions (Ilbert et al. 2009) and apply the cuts on size, 
magnitude and redshift summarized in Tab. 1 . The size cut ap- 
proximately corresponds to a measured size (before PSF de- 
convolution) of 1.4 times the PSF size. Note that in the shear 
analysis, a significantly stricter cut {d > 3.6) is employed. 
The magnitude cut is mostly determined by the reliability of 
the photometric redshifts. Finally, the redshift cut excludes 
a small portion of the sample which does not contribute sig- 
nificantly to the lensing signal-to-noise. We obtain a sample 
of A^src = 250,500 galaxies with a median redshift of 1.22. In 
the lensing analysis, we also impose the conditions z, _ > zl, 
Zi > zl + 0.2 for each source-lens pair. 

We apply the fitting procedure described in the last section, 
using overlapping redshift slices of half- width Az = 0.25 sep- 
arated by Az. Fig. 1 illustrates the measured x(z), m{z) as well 
as (Ti.int, (Tm.int- We fix /9 = -0.6 to improvc the stability of the 
fit. Allowing p to vary does not significantly improve the fit. 
Fig. 2 shows Ex and Note that < 0: the number of faint, 
small galaxies passing the cuts for k >Q reduces the average 
size, counteracting the lensing effect on the sizes themselves. 
The analogous holds for magnitudes, where £,„ > 0. This ef- 
fect can in principle be avoided by cutting on a quantity not 
affected by lensing (e.g., surface brightness) instead of mag- 
nitude and size. We leave the exploration of this for future 
work. 

While the estimators in Sec. 2 were designed to rely as 
much as possible on data alone, we do need simulations to de- 
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Fig. 2. — From top to bottom: a) Flux lensing efficiency q an a function 
of magnitude m. b) Lensing bias correction e, for galaxy (log) sizes. The 
effective size lensing efficiency is given by fj = rj + Sx- c) Lensing bias effect 
£,„ for magnitudes. Correspondingly, q = q + e,,, . 

termine the lensing efficiencies rj, q. To simulate ACS data, 
we use the s image software package (Massey et al. 2004; 
Ferry et al. 2008; Dobke et al. 2010), as previously adopted 
for the Shear TEsting Programme (STEP; Massey et al. 2007) 
and to calibrate the COSMOS shear analysis (Leauthaud et al. 
2007). The parameters were tuned to mimic the galaxy source 
counts and correlated background noise properties of the re- 
duced and stacked COSMOS images. We model the intrin- 
sic morphologies of source galaxies via shapelets (Refregier 
2003a; Massey & Refregier 2005). We have adapted s image 
to enable the superposition of an input convergence n^im and 
produced images at 15 values of Ksim G [-0.2,0.3]. 

To determine how each simulated galaxy responds to mag- 
nification, we match galaxies in the lensed and unlensed 
(usim = 0) simulations. We measure the average x and m as 
a function of the input Ksim = -0.1 . . .0.1. Then, we estimate 
?7 through 

11= , (25) 

K 

and correspondingly for q. We have found that for galaxies 
passing cuts, 7] is within [0.9,1.1] and consistent with 1 for 
the whole sample (recall that we use a PSF-deconvolved size 
estimator). 

For the lens sample, we use X-ray selected groups since 
these have been well studied using shear (Leauthaud et al. 
2010). The entire COSMOS region has been mapped by 
XMM-Newton, while the central region is covered by Chan- 
dra observations (Hasinger et al. 2007; Cappelluti et al. 2009; 
Elvis et al. 2009). These X-ray data have been used to con- 
struct a group catalog containing 211 extended X-ray sources, 
spanning the range < Zgr < 1 (Finoguenov et al. 2007; Leau- 
thaud et al. 2010). Group membership is assigned to galaxies 
in the COSMOS catalog based on the photometric redshifts 

TABLE 1 

Summary of cuts applied to the source galaxy sample. 



m 


< 25.8 mag 
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> 2 pixels 
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Fig. 3. — Mean reconstructed convergence k in simulations with magnifi- 
cation, as a function of tlie input Ksi,,,. The thin solid line shows the identity 
while the thick solid line shows a quadratic fit. 

and the red sequence. We include only groups with reliable 
optical associations which have more than 3 members, and we 
exclude close neighboring systems and groups near the edges 
of the field or masked regions (FLAG_INCLUDE=1 in George 
et al. (2011a)). We use the location of the most massive group 
galaxy (MMGG) located within the NFW scale radius of the 
X-ray centroid as the group center (MMGGscaie as defined in 
George et al. (201 lb)). George et al. (201 lb) show that this is 
likely the most reliable estimate for the group centers. We re- 
strict the redshift range of groups to 0.2-0.6, close to the peak 
of the lensing sensitivity of the source sample. This selec- 
tion also reduces the potential impact of catastrophic redshift 
errors which become more prevalent at > 1.5, and yields 
A^gr = 61 groups. 

4. RESULTS 

Before applying our estimators to data, we performed the 
calibration described in Sec. 2 on the unlensed simulations, 
and applied the estimators to the lensed simulations as a con- 
sistency check. Since the simulations do not contain redshift 
information, we set F(zi) = 1. Fig. 3 shows the mean recon- 
structed K as a function of the true Ksim, along with a quadratic 
fit, 

(Krec) =P0 + PlKsim + P2Ksim- (26) 

Note that we expect some quadratic correction in this rela- 
tion, since our estimators neglect all second- and higher-order 
terms. We found pi = 0.95 ±0.05, indicating that our estima- 
tor is close to unbiased, and p2 ~ 0.6. 

Turning to data, we collect source galaxies in bins of physi- 
cal transverse radius r around the A^gr X-ray groups, and apply 
the estimator for S. The second order corrections amount to 
at most ~' 10% at the smallest radii, as shown in Leauthaud 
et al. (2010), and can be neglected given the errors of this 
measurement. Since x(z), m{z) are measured in finite redshift 
bins without any redshift weighting while we apply redshift 
weighting in the E measurement, we need to subtract a con- 
stant offset Soff from the measurement. Soff is measured by 
pairing all galaxies passing cuts in the COSMOS field with 
each group, yielding Soff = 151.9Mo/pc-. This is equiva- 
lent to estimating (S(/"))gi-oups - (S)cosmos, where both aver- 
ages employ the same redshift weighting. Note that a circle 
of r 1 Mpc drawn around each group covers a significant 
fraction of the COSMOS field; thus, we necessarily subtract 
some of the signal in this way. If we instead restrict the source 
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Fig. 4. — Projected surface mass density S (filled circles) measured around 
groups using the estimator described in Sec. 2. The open circles show AS 
measured using shear for the same group sample. The solid and dotted 
lines show E and AE, respectively, for the joint best-fitting NFW model 
(lg(M/MQ) = 13.43; r < 1 Mpc). Note that associated large-scale structure 
contributes significantly to S(r) for r > 1 Mpc. 

sample used for the Eoff estimate to galaxies separated by at 
least 1 Mpc in physical transverse distance from each group 
in the sample, Sotf decreases by 4.4Mq/pc-, thus increasing 
Ti{r) uniformly by this amount. 

The measured I](r) is shown in Fig. 4. The error bars are 
those returned by the Gaussian estimators. Measurements of 
the variance of S applied to random sets of source galaxies 
indicate that these errors are correct to within ^15%. For 
r < 1 Mpc, the errors for different radial bins are independent. 

We fit an NFW profile (Navarro et al. 1997) with fixed con- 
centration c = 4 to the measurement, calculating I](r) as de- 
scribed e.g. in Schmidt & Rozo (2011). Furthermore, we re- 
strict the fitting to r < 1 Mpc, beyond which the contribution 
from associated large-scale structure becomes important (see 
below). This yields a best-fit mass of lg(M/M0) = 13.25 ± 
0.16 (statistical error), corresponding to a detection signifi- 
cance of ~ 3.8(7. This corresponds to ^ 38% of the signal- 
to-noise of shear Note that the precise value also depends on 
the definition of the offset Soff. The best-fit mass is consis- 
tent with the estimate from shear, lg(M/M0) = 13.49 ±0.07 
(Leauthaud et al. (2010); again for r < 1 Mpc). 

In principle, a measurement of E in addition to AS can 
break the mass-concentration degeneracy, since the two quan- 
tities have different radial dependencies (Rozo & Schmidt 
2010). Due to the limited signal-to-noise of this first measure- 
ment, a joint fit to the S measurement together with AS from 
shear only yields an improvement of ~' 10% in the area of the 
error ellipse in the mass-concentration plane. The S measure- 
ments at r > 1 Mpc are consistent with an estimate of the asso- 
ciated large-scale structure contribution (two-halo term) using 
the halo model, which yields S2/, ~ 10Mq/pc^ with weak r- 
dependence for r < 5 Mpc. The corresponding contribution 
to the shear is much smaller, AS2/, < 1 Mo/pc^, further ex- 
emplifying the complementarity of the two measurements. A 
detailed modeling of these data will be the subject of forth- 
coming work. 

The main systematic uncertainty of the S measurement is 
due to the lensing bias effect. We have quantified the impact 
of the measurement errors on e^, e„, by varying nicm, x^at, re- 
sulting in an overall multipHcative systematic error on S of 
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20%. Residual effects from PSF and background variations 
have been found to be negligible. Increasing the minimum 
size cut from 2 to 3 pixels yields compatible results, as do 
estimators using sizes and magnitudes separately. 

The systematic uncertainty due to photometric redshift er- 
rors is similar to that of the shear (Leauthaud et al. (2010); 
note that we have employed the same source-lens separation 
criteria as in that paper). Given the lower signal-to-noise of 
the magnification measurement, this systematic is negligible. 

5. SUMMARY AND DISCUSSION 

We have presented the first measurement of weak lensing 
magnification using galaxy sizes and fluxes. Measurements 
of galaxy sizes and magnitudes come for free with weak lens- 
ing surveys. We find a signal-to-noise only a factor of 2 - 3 
less than shear; thus it is worth pursuing this kind of indepen- 
dent weak lensing measurement. Magnification constitutes a 
different observable than shear, and the non-local relationship 
between the two can be used to break degeneracies e.g. when 
measuring halo profiles. 

Furthermore, the systematics affecting this measurement 
are largely independent from those of the shear: to first or- 
der, sizes are only sensitive to the size of the PSF, while mag- 
nitudes are PSF-independent. On the other hand, a knowl- 
edge of the shape of the PSF is necessary for shear measure- 
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ments. This allows us to use less restrictive cuts. Conversely, 
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sary for the magnification measurement, while this effect is 
not present in the shear at linear order (Schmidt et al. 2009a). 
Nevertheless, given sufficiently detailed simulations, this is a 
solvable problem. 

While the measurement presented here is limited by signal- 
to-noise and statistics, it is clear that for future, much larger 
surveys, robust size and flux estimators are essential to a suc- 
cessful application of this method. Given the already promis- 
ing results, we believe a dedicated effort to develop such es- 
timators for weak lensing will be worthwhile. Finally, our 
results show that there is significant lensing information in 
photometric galaxy samples in addition to the shear alone; 
this should be taken into consideration in the design of future 
weak lensing experiments. 
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